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ABSTRACT 

We investigate the correlation between the mass of the supermassive black holes 
(SMBHs) and metal abundance, using existing data sets. The SMBH mass Mhh is 
well correlated with integrated stellar feature of Mgb. For 28 galaxies, the best-fit 
Mb/j-Mgb relation has a small scatter, which is an equivalent level with other well- 
known relation, such as a correlation between the stellar velocity dispersion and the 
mass. An averaged iron index (Fe) also positively correlates with Mhh, but the best- 
fit Mbh-(Fe) relation has a larger scatter. The difference comes from the synthesis 
and evolution mechanisms, and may be important for the SMBH and star formation 
history in the host spheroid. 
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1 INTRODUCTION 

Recent observations of nearby massive spheroids (ellipticals, lenticular and spiral bulges) have established that supermas- 
sive black holes (SMBHs) are present in the nuclei of galaxies. The mass Mbh of SMBHs ranges from 10 6 to 1O 9 M . 
Several statistical correlations with characteristic parameters of the host spheroids are explored: relatio n with the hum 



nosity L{ KOTmend y feRichstone! 1995 ; Marconi fe Hunt 2003 



Graham 



2007) , the mass M s of the sph eroids) Magomanetal 



the light concentration or Sersic index n ([Graham et al.l l200ll ; iGraham fe Driver! 120071 ) . gravitational binding energy E, 



199^ : iHaring fe RbJl2004l ). the stellar veloci ty dispersion a [ Gebhardt et al.l 2000l: iFerrarese et alJbOQCl : iTremaine et al.ll2002h 



Aller fe Richstone 2007 ) and so on. The recent updated survey for various galaxy parameters is given in Aller fe Richstone] 
2007] ). If the fitting models are good, the correlations may provide some clues to the coevolution of S MBHs and the hos t 



spheroids. Some tentative theories have been proposed to clarify the origin of these correlations, e.g., Silk fe Rees (1998) 



Umemural (2001) 



The galaxies have some aspects as an assemble of stars, gases and dark matters. The quantities M s , a, E g represent 
the dynamical aspect, while L, n represent the photometric one, although they are closely related each other. Very little 
attention is paid to the chemical one. The observation is not so easy and therefore the correlation between SMBH mass and 
spheroid chemical property is not well studied so far. Both the metal abundance and black hole mass increase with time and 



chemical parameter u sing Lick/IDS absorption line indices, e 



studied in literatures (|Trager et al.ll200ol : iDenicolo et al.l 1200 



.. Worthev et al. di 


994), 


Kuntschner et al.1 


2006) 



SMBH is reliable. Using the published data, we investigate the correlation with SMBH mass in this paper. 

As a sample we select only 28 galaxies, for which SMBH masses are measured by dynamical methods. In section 2, we 
describe the sample of galaxies and criterion of our choice. In section 3, we derive the correlation between Mbh and a, and 
that between Mbh and the B-band magnitude Mb from our sample. These correlations are known to be good, and are used in 
order to check no serious bias inv o lved in the sam ple. Our results for Mbh-& and Mbh-Ms relations are not so different from 
those by Aller fe Richstone] ( 2007 ); Graham 1 2007 ). although the source sample is not the same. We also show the correlation 
between black hole mass and metal abundance index using our sample. Finally, a discussion of our results is given in section 
4. 
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2 THE SAMPLE 

TABLE 1 Galaxy Parameter Used in The Fits 

Galaxy Type" (low, high) Reference a b Mb c Mgb (Fe) Reference 

(Mq) (km s~ x ) (mag) (mag) (mag) 



(1) 


(2) 




(3) 


(4) 


(5) 


(6) 






(7) 






(8) 




(9) 


NGC 221 


SO 


2.5E6 


(2.0,3.0) 


1 


75 


-15 


3 





099 


± 





001 





078 


± 





001 


19,20 


NGC 224 


Sb 


1.4E8 


(1.1,2.3) 


2 


160 


-19 








175 


± 





002 





087 


± 





001 


20 


NGC 821 


El 


8.5E7 


(5.0,12) 


3 


209 


-20 


4 





159 


± 





002 





086 


± 





001 


19,20,21 


NGC 1023 


SBO 


4.4E7 


(3.9,4.9) 


4 


205 


-18 


4 





162 


± 





003 












21 


NGC 2974 


El 


1.7E8 


(1.5,1.9) 


5,6 


233 


-20 


9 





159 


± 





002 





084 


± 





002 


19,21 


NGC 3115 


SO 


1.0E9 


(0.4,2.0) 


7,8 


230 


-20 


2 





177 


± 





002 


o 


096 


± 





001 


19,22 


NGC 3245 


SO 


2.1E8 


(1.6,2.6) 


9 


205 


-19 


6 





161 


± 





006 





077 


± 





002 


19 


NGC 3377 


E5 


1.0E8 


(0.9,1.9) 


10 


145 


-19 








142 


± 





002 





076 


± 





001 


19,20,21 


NGC 3379 


El 


1.4E8 


(0.6,2.0) 


11 


206 


-19 


9 





168 


± 





001 





084 


± 





001 


19,20,21 


NGC 3384 


SBO 


1.6E7 


(1.4,1.7) 


10 


1 13 


-19 








154 


± 





002 





075 


± 





001 


19,21 


NGC 3414 


SO 


2.5E8 


(2.2,2.8) 


5,6 


205 


-20 








170 


± 





003 


o 


081 


± 





001 


19,21 


NGC 3608 


E2 


1.9E8 


(1.3,2.9) 


10 


182 


-19 


9 





162 


± 





002 





085 


± 





001 


19,20,21 


NGC 4261 


E2 


5.2E8 


(4.1,6.2) 


12 


315 


-21 


1 





186 


± 





002 





088 


± 





001 


19,20 


NGC 4374 


El 


4.6E8 


(2.8,8.1) 


13,14 


296 


-21 


4 





168 


± 





001 





082 


± 





002 


19,21,22 


NGC 4459 


SO 


7.0E7 


(5.7,8.3) 


15 


186 


-19 


1 





139 


± 





003 












21 


NGC 4473 


E5 


1.1E8 


(0.3,1.5) 


10 


190 


-19 


9 





168 


± 





002 





088 


± 





002 


21,22 


NGC 4486 


E0 


3.0E9 


(2.0,4.0) 


7,16,17 


375 


-21 


5 





199 


± 





002 





086 


± 





002 


21,22 


NGC 4552 


E 


5.0E8 


(4.5,5.5) 


5,6 


252 


-19 


2 





185 


± 





002 





084 


± 





001 


20,21 


NGC 4564 


SO 


5.6E7 


(4.8,5.9) 


10 


162 


-17 


4 





173 


± 





003 












21 


NGC 4596 


SBO 


7.8E7 


(4.5,11.6) 


15 


152 


-20 


6 





150 


± 





008 












23 


NGC 4621 


E5 


4.0E8 


(3.6,4.4) 


5,6 


211 


-19 


5 





178 


± 





002 





088 


± 





002 


21,22 


NGC 4649 


El 


2.0E9 


(1.4,2.4) 


10 


385 


-21 


3 





194 


± 





002 





085 


± 





001 


20 


NGC 4697 


El 


1.7E8 


(1.6,1.9) 


10 


177 


-20 


2 





146 


± 





002 





078 


± 





001 


20 


NGC 5813 


El 


7.0E8 


(6.3,7.7) 


5,6 


230 


-20 


9 





167 


± 





001 





078 


± 





001 


19,20,21 


NGC 5845 


E3 


2.4E8 


(1.0,2.8) 


10 


234 


-18 


7 





169 


± 





003 





093 


± 





002 


19,21 


NGC 5846 


E0 


1.1E9 


(1.0,1.2) 


5,6 


238 


-21 


2 





179 


± 





001 





085 


± 





001 


19,20,21,22 


NGC 7052 


E4 


3.3E8 


(2.0,5.6) 


7,18 


266 


-21 


2 





182 


± 





002 





080 


± 





002 


20 


NGC 7457 


SO 


3.5E6 


(2.1,4.6) 


10 


67 


-17 








101 


± 





003 












21 



NOTES.-Col.(l):Galaxy name. Col.(2):Galaxy type. Col.(3):SMBH mass. Col.(4):Reference for col.(3). Col. (5) Effective stellar bulge 
velocity dispersion. Col. (6): Absolute bulge B-band magnitude. Col.(7),(8)Line strength index measurements of the luminosity weighted 
spectrum within r e /8 for the Mgb and (Fe) indices. Col.(9):Reference for col. (7), (8). 

a Morphological type from NED e xcept for NGC 224 and NGC 4564. The type given here for those galaxies is from [Graham] ^2002) 
for NGC 221, and from lTruiillo et alj j2004l) for NG C 4564. 

b The value is taken fromlTremaine et"afl (120021 1 except for NGC 2974, NGC 3414, NGC 4552, NGC 4621, NGC 5813, NGC5846lHul 
1120081) , NGC 4374 iMarconi fe Huntl j2003h. 

c Th e value is taken fromlGrahaml ||2007') or extracted from IMarconi fe Hum] J2003T ) and the Third Reference Catalogue of Bright 
Galaxies Ide Vaucouleurs et al.l dl99lT) . 

REFERENCES,- (T lVerolme et al.1 [120021 ) f2 feender et alj j2005t) f3 lRichstone et al.l ||2004T ) f4 lBower et alj feOOlft (5 lCappellari et al.l 
i2007h (6 ~lHulj2008h (7 lTremaine et al.l j2002h (8 lKormendv et al.Nl996l) Q lBarth et alJfeOOlD ( lO lGebhardt et al.l (120031) (Tl lShapiro et alj 
J2006I) (12 lFerrarese. Ford fe Jaffd jl996h (T3 lMacieiewski fe Binnevl J200ll) fU lKormendv et alJfeOOll) fl5 jSarzi et aLlfeOOll) (16 lHarms et alj 
Jl994l) (17 lMacchetto et al.Nl997l ) fl8 jvan der Marel fe van den Boschl dl99St) fl9 lDenicol6 et al.l j2005l) (20 lTrager et al.ld2000l) (2l lKuntschner et alj 
l|2006l) f22 lHowelll j2005l ) f23 jPeletier et alj d2007tl 



Our sample consists of 28 nearby galaxies, which are listed in Table 1. We put two criteria to pick up the sample from 
literatures. First, the BH masses are measured by good spatial resolution. The sources are limited to Mb h > a 2 r res /(2G) 



where a is stellar velocity dispersion, r res the instrumental spatial resolution. This condition is used by e.g., Marconi & Hunt 



2005 ); Hu (200.8). Second, metal indices of Mg or Fe are measured. Most of our sample corresponds to that of Tremaine et al 



20021 ). 21 galaxies are drawn from their table and 7 galaxies a re added: NGC 2974, NGC 3414, NGC 4552 , NGC 4621, 



NGC 5813 and NGC5846 from ICappellari et al.1 (|2007l ); IHuI (120081). and NGC 4374 from iMacieiewski fe Binnevl |200ll) . The 
mass es of SMBHs are u pdated for 3 galaxies, NGC 224 j|Bender et al.l I2005T I. NGC 821 (|Richstone et al.1 booi ) and NGC 
3379 (|Shapiro et al.1 I2OO6I ). The galaxy morphology class is also liste d in Table 1 from t he NASA Extr agalactic Data base 
(NED). However, some corrections are added to the data: NGC 4564 (|Truiillo et all Exil ) and NGC 221 (|Grahamll2002l ) are 
now recognized as SO ga laxies. In Table 1 , effec tive velocity dispersion a is tabulated as an indicator of the dynamical parameter 
of spheroids. Following iTremaine et alj (|2002l ). the relative errors in a are estimated as 5 % (0.021 dex). As the photometric 
parameter of spheroids, the B-band magnitude Mb is tabulated. Following iGrahaml l|2007l ). the relative errors in Mb are 0.3 
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mag. For t he chemical paramete r of spheroids , we u se spectral absorption line indices in the visual as defined by the Lick 
group, e.g., iBurstein et all (|l984l ); iFaber et all 1119851). Lick indices ha ve proven to be a useful tool for the derivation of ages 
and metallicity of unresolved stellar populations (|Denicol6 et al,ll2005l ). In this paper, we use Mgb, Fe5270, Fe5335 indices, 
measuring respectively the strength of Mgl at A ~ 5156 - 5197A, Fel at A ~ 5246 - 5286Aand Fel at A ~ 5312 - 5352A. These 
indices in mag are calculated by the standard equations: 



EW„ 



-2.5 loe 



J^[F(X)/C(\)]d\ 



A2 — Ai 



(1) 



\ r — A 



A- A„ 



(2) 



where 

A r A^ At- A^ 

and Ai, and A r are the mean wavelength in the blue and red pseudo-continuum intervals, respecti vely. We have adopte d the 
spectral pseudocontinua and band passes of the Mgb, Fe5270, Fe5335 Lick/IDS indices defined in lWorthev et al.l ([19941 ). We 
use a combined "iron" index, (Fe) defined by 

(Fe) = i(Fe5270 + Fe5335). (3) 

This index has smaller random error than either that of Fe5270 or that of Fe5335 (|Gorgas. Efstathiou fe Aragon-Salamanca 
1990). These index values used from the literatures, are the measurements within the central region r < r e /8, where r e is 
effective radius of the spheroids. The properties of the spheroids and SMBHs are listed with their references in Table 1. 



3 RESULTS 



In order to fit the data to the linear relation y = a + bx, we use a version of the routine FITEXY IjPress et al.lll992l ) modified 
by iTremaine et al.l ( 20021 ) . iNovak et alj ( 20061 ) showed that this algorithm is the most efficient and least biased among a set 
of algorithms explored by them. The best-fit values a and b are calculated from data (xi,yi)(i = 1, • • • , N) by minimizing the 
quantity 



N 

2 \ " 

X = 



bxi 



^ 5y1 + e 2 + b 2 8x 2 t ' 

i — 1 



(4) 



where 5xi and 5yi are the measurement uncertainty in Xi and yi. We use a constant value of 5 % (0.021 dex) for velocity 
dispersion a, and 0.3 mag for the B-band magnitude Mb- The uncertainty for the metal abundance is listed in Table 1, and is 
used for each source. The intrinsic scatter e in the linear relation is calculated when the reduced chi-squared value X 1 /(N — 2) 
is equal to 1, after many trails to fit it for various e. The upper and lower uncertainties of e are estimated by the values for 
X 2 /(N - 2) = 1 ± J2/N. 



3.1 Dynamical and Photometric Properties 



In order to check the dependence on the sampling data, we derive M^-u and M^-Mb relations from our sample. First, for 
the dynamical aspect, we show a relation between the SMBH mass and the stellar velocity dispersion derived from Table 1. 
Our result of M^h-o relation is given by 



log(M bh /M Q ) = (3.98 ± 0.35) log(a/200) + 8.23 ± 0.06 



(5) 



with e = 0.25ln'f!|dex in log Mbh- This relation is consistent with previous relations. For example, recent updated one by 



Aller fc Richstond (|2007D using a sample of 23 galaxies is 



log(M bh /M Q ) = (3.79 ± 0.32) log(cr/200) + 8.16 ± 0.06 

with e — 0.23 dex. Figure 1 shows our Mhh-o relation ((5JI in comparison with (BJ. 
Next, for the photometric aspect, we derive the M^h-Mb relation. Our result is 

\og(M bh /M Q ) = (-0.42 ± 0.05)(A/ S + 19.5) + 8.14 ± 0.08 



with e = O^e+g^deaH lGrahaml |20071 ) obtained 
\og{M bh /M Q ) = (-0.40 ± 0.05) (M B + 19.5) + 8.27 ± O.i 



(6) 

(7) 
(8) 



1 We have also derived the relation of with the K-band magnitude, but there is no significant difference in the scatters of both 
relations. 
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log a 

Figure 1. Relationship between SMBH mass and effective 
stellar velocity dispersion for 28 galaxy samples. The solid 
line represen ts our best fit relation iQ. For a comparison, the 
relation by (All er fc Richstonel 120071) is shown by a dashed 
line. Elliptical galaxies are denoted by open circles, lenticulars 
and spirals by filled circles. 



10 




g I i i i i i i I 

-15 -16 -17 -18 -19 -20 -21 -22 
Mo 



Figure 2. Relationship between SMBH mass and absolute 
B-band luminosity of the spheroid for 28 galaxy samples. The 
solid line represents our best fit relatio n Q. For a comparison, 
the relation found by l lGrahamll2007n is shown by a dashed 
line. The symbols are the same as Figure 1. 



i 

CD 



o 




0.1 0.12 0.14 0.16 0.18 0.2 0.07 0.075 0.08 0.085 0.09 0.095 0.1 

Mgb <Fe> 

Figure 3. Relationship between SMBH mass and Mgb index Figure 4. Relationship between SMBH mass and (Fe) index 
value within a circular aperture of r e /8 for our samples. The within r e /8 for our samples. The symbols and line are the 
symbols and line are the same as in Figure 1. same as in Figure 1. 

with e = O.SOlg Qgdex. Although intrinsic scatter in our relation is a little high, both two relations agree within the uncer- 
tainties. Both relations (0 and |[8j are plotted in Figure 2. By comparing with previous results, we may say that there is no 
serious bias in our sampling data. 



3.2 Chemical Property 

In Figure 3, we plot M&j, as a function of the integrated stellar feature of Mgb for the 28 galaxies. It is found that the mass 
M\,h significantly correlates with the index Mgb. The best fit relation is obtained as 

\og{M bh /M Q ) = (27.23 ± 3.00)(Mg6 - 0.16) + 8.13 ± 0.07 (9) 

with e — 0.32lg Q|dex. The scatter is not so bad as that of well-fitted Mj^-cr and M^^-Mb relations. The index Mgb is a good 
indicator, but other indicator of the metal abundance is not so good. We have tested the correlation with different index. For 
instance, we show M b h as a function of (Fe) index in Figure 4. The best fit relation is 

\og{M bh /M Q ) = (60.25 ± 26.01) ((Fe) - 0.083) + 8.34 ± 0.13 (10) 

with e = 0.58^0 07 dex. This relation has much larger scatter than that of eqs.(5)-(9). It is clear from Figure 4 that the index 
(Fe) is not good one. Summarizing the relationships obtained in our own study, we find that M^-cr is the best, M^-Mgb and 
Mth-Ms are moderate, and Mi,h-(Fe) is the worst. 



4 DISCUSSION 

We have shown that a good correlation between SMBH mass and Mgb index value. The best-fitting M&ft-Mgb relation has 
small intrinsic scatter 0.32 dex which is comparable one in other strong correlations found so far. Such a new correlation in 
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the chemical aspect is expected through other relations in dynami cal/ ph o tomet ric aspects. The metal abundance is roughly 
correlated with total stellar mass, absolute magnitudes etc. e.g., IPagell dl997h. In part i cular, there is a re markably tight 



relation between Mg2 and the central velocity dispersion of stars (|Terlevich et alj Il98ll ; iBender et al.lll993h . The positive 
correlations have been found between some quantities characterizing the hosts and SMBH mass. A positive correlation is 
also suggested bet ween the S MBH ma ss and the met allicity derived from emission line ratios in 578 AGNs spanning a wide 
range in redshiftsl Warner. Hamann fc Diet rich 2003). Our result of Mgb for nearby galaxies is more tight. Thus, positive 
correlation is expected, but the degree was not clear beforehand. It was not clear which indicators of the metal abundance 
strongly correlate with the SMBH mass. The index (Fe) correlates with it, but the intrinsic scatter is not so small as that of 
Mgb. Heavy elements Mg and Fe are synthesized by two different types of supernovae (Type la and II), with different time 
scales. This evolutionary difference causes the scatter in the correlations. 

The growth of black hole mass is mainly determined by the accretion rate and the lifetime of the activity. The environ- 
mental factors near the central region of galaxies may partially be affected by some global quantities, such as the mass and 
size of the host. If the mass Mf,h of SMBH is determined solely by the spheroid mass M s as Mbh = eM s (e ~ 10 -5 ), then we 
have Mbh °^ M s oc L oc a 4 , where we assume that a constant mass-to-light ratio and the Faber- Jackson relation in elliptical 
galaxies hold in the spheroids. The features in the host spheroids are transferred to the relations with the black hole mass, 
Mbh-L an d Mbh-cr relations. Ot her features in the hosts, binding energy, light concentration and so on also give some relations 
with Mbh- Bender et al. ( 1993h discuss that the strength of Mg2 is determined not only by the global mass M s , but also by 
the local stellar density, which is related with the star-formation rate etc. It is therefore important to examine the correlation 
of A4bh with the other physical quantities of the spheroid, in addition to M 3 . The metal abundance is a tracer for integrated 
stellar populations. The tight relation Mbh-Mgb, which is discovered here but is still tentative, may be useful for the better 
understanding of the coevolution of SMBH and host spheroid, if it is not accidental. Further work is needed to clarify whether 
this relation is fundamental or not. 
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